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l.  Introduction 

We  have  developed  a  bulk  PBL  model  with  a  simple 
internal  vertical  structure  and  a  simple  second-order  closure, 
designed  for  use  as  a  PBL  parameterization  in  a  large-scale 
model.  The  turbulent  fluxes  in  the  interior  of  the  PBL  are 
parameterized  using  the  "convective  mass  flux"  concept 
introduced  by  Arakawa  (1969)  and  adopted  in  many 
subsequent  studies.  Fig.  1  summarizes  the  conceptual 
pedigree  of  this  second-order  bulk  model"  (S.O.B.),  relative 
to  earlier  models  used  in  boundary  layer  and  cumulus 
parametenzations.  A  sketch  of  the  model  s  formulation  is  as 
follows: 
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Figure  1  Diagram  summarizing  the  relationship  of  the 
present  model  to  earlier  models  used  in 
boundary-layer  and  cumulus  parametenzations. 


The  basic  framework  of  the  model  is  shown  in  Fig.  2. 
The  level  just  above  the  PBL  top  is  denoted  by  subscript  B~. 
while  the  Earth's  surface  is  denoted  by  S-.  We  define  an 
infinitesimal  "entrainment  layer"  just  below  the  PBL  top.  and 
an  infinitesimal  "ventilation  layer"  just  above  the  Earth  s 
surface.  The  concept  of  an  entrainment  layer  is  motivated  by 
the  observations  of  Caughey  er  al.  ( 1982)  and  Nicholls  and 
Turton  (1986),  who  described  it  as  a  thin  region  of  weak 
organized  vertical  motions  and  vigorous  small-scale  mixing. 
The  ventilation  layer  is  more  conventionally  known  as  the 
surface  layer. 

B+ _ 


Figure  2:  Diagram  illustrating  the  assumed  structure  of  the 

PBL.  The  interior,  which  is  represented  by  two 
layers,  is  bounded  above  by  a  thin  entrainment 
layer  and  below  by  a  thin  ventilation  layer 
Convective  circulations  occur,  with  ris.r. z 
branches  occupying  fractional  area  (7  The 
ascending  and  descending  branches  have 
different  thermodynamic  soundings  jnd. 
therefore,  different  cloud  base  levels. 


The  depth  of  the  PBL  (in  terms  of  pressure)  is 
prognosttcallv  determined,  and  is  denoted  by  <5p«.  The 
generic  variable  w  represents  a  prognostic  intensive  scalar 
such  as  the  dry  static  energy,  the  the  mixing  ratio  of  water,  or 
a  component  of  the  horizontal  wind.  Area-averaged  values  of 

ip  are  denoted  by  yf .  The  turbulent  flux  of  ip,  denoted  by  Fy, 
is  defined  at  the  top  of  the  surface  layer,  and  at  a  level  just 
below  the  PBL  top.  These  levels  are  denoted  by  subscripts  S, 
and  B.  respectively.  The  turbulent  momentum  flux  is  also 
defined  levels  S  and  B. 

An  entrainment  mass  flux.  E.  carries  mass  across  the 
PBL  top.  and  is  closely  related  to  the  turbulent  fluxes  near  the 
PBL  top.  Correspondingly,  a  ventilation  mass  flux,  V,  is 
associated  with  the  surface  fluxes:  in  conventional  parlance,  V 
is  the  product  of  the  surface  air  density,  the  surface  wind 
speed,  and  a  transfer  coefficient.  For  both  the  entrainment  and 
ventilation  layers,  the  model  incorporates  diagnostic  balances 

for  mass.  ip.  and  ip'2.  The  entrainment  and  ventilation  layers 
are  assumed  to  be  thin  enough  so  that  such  balance  conditions 
are  appropriate.  Within  the  ventilation  layer,  the  turbulent 
fluxes  have  to  be  earned  by  small  eddies,  since  the  organized 
vertical  motions  associated  with  the  convective  circulations 
must  vanish  there.  The  ventilation  layer  is  assumed  to  be  thin, 
in  the  sense  that  the  turbulent  fluxes  at  its  top  are 
approximately  equal  to  those  at  the  surface.  Similarly,  within 
the  entrainment  layer  the  organized  vertical  motions  associated 
with  the  convective  circulations  become  negligible,  so  that 
smaller  eddies  must  again  carry  the  turbulent  fluxes.  The 
entrainment  layer  is  assumed  to  be  thin  in  the  sense  that  the 
turbulent  fluxes  at  its  base  are  approximately  equal  to  those  at 
the  PBL  top. 

The  vertically  integrated  TKE  of  the  PBL  is 
prognostically  determined,  and  is  denoted  by  ey.  In  addition, 

the  vertically  integrated  scalar  variances,  (ipTw,  are 
diagnostically  determined.  Because  the  model  makes  use  of 
these  second-order  turbulence  variables,  it  can  be  viewed  as  a 
highly  simplified  second-order  closure  model  with  very  coarse 
vertical  resolution.  The  coarse  vertical  resolution  is  made 
palatable  (and  feasible)  by  the  use  of  an  explicit,  prognostic 
PBL  depth  The  fine  vertical  structures  that  typically  occur 
near  the  PBL  top  and  the  Earth  s  surface  are  parametrically 
included  in  the  model,  by  methods  similar  to  those  used  in 
mixed- layer  models.  This  parametric  representation  of  the 
vertical  structures  of  the  entrainment  and  ventilation  layers  is 
an  alternative  to  explicitly  resolving  these  thin  layers;  the  latter 
approach  is  followed  in  conventional  higher-order  closure 
models  and  is  the  primary  reason  that  such  models  require 
high  vertical  resolution. 


By  assuming  that  the  ventilation  and  entrainment  layers 
at  the  lower  and  upper  edges  of  the  PBL  are  dominated  by 
small-scale  turbulence  in  quasiequilibnum.  boundary 
conditions  are  developed  for  the  rising  and  sinking  branches  of 
the  convective  circulations,  and  also  for  the  scalar  variances 
associated  with  the  convective  circulations.  The  convective 
mass  flux  and  the  fractional  area  covered  by  updrafts  are 
diagnosed  by  the  model.  Fractional  cloudiness  occurs  when 
the  ascending  branches  are  saturated  and  the  descending 
branches  are  not. 

The  convective  mass  flux  is  determined  using  the 
predicted  TKE.  In  this  same  spirit,  we  determine  the  surface 
fluxes  using  a  modified  bulk  formu'a  in  which  the  square  root 
of  the  TKE  takes  the  place  of  the  wind  speed.  This  allows  a 
natural  representation  of  free  convection,  and  also  minimizes 
the  stability  dependence  of  the  effective  transfer  coefficient. 
Finally,  following  Siems  et  at  (1990).  the  entrainment  rate  is 
also  assumed  to  be  proportional  to  the  square  root  of  the  TKE; 
ihe  proportionality  factor  depends  on  the  inversion  Richardson 


number,  and  also  on  an  additional  parameter  that  represents  the 
effects  of  evaporative  cooling  when  clouds  are  present. 


A  detailed  description  of  the  model's  formulation  is 
given  by  Randall  and  Shao  1 1990).  Here  we  describe  only  the 
method  used  to  determine  the  fracnonal  area  covered  by  risina 
motion. 


2.  A  method  to  determine  cr 

A  key  to  determining  the  fractional  cloudiness  is 
finding  the  value  of  cr.  the  fractional  area  covered  bv  nsinc 
motion.  Previous  applications  of  the  mass  flux  model  to 
boundary  layer  clouds  have  not  included  a  method  to 
determine  cr.  A  brief  summary  of  our  approach  to  this 
problem  is  as  follows: 

W'e  assume  that,  in  the  interior  of  the  PBL.  the 
turbulent  fiuxes  are  entirely  due  to  the  convective  circulations, 
with  rising  branches  covering  fractional  area  cr.  and  sinking 
branches  covering  fractional  area  /  -  a  Observations  based 
on  conditional  sampling  methods  suggest  that  O’ is  typically 
less  than  1/2  for  the  clear  convective  PBL  (e.g..  Lenschow 
and  Stephens,  1980),  and  greater  than  1/2  for  the  cloud-topped 
PBL  (e.g..  Nicholls,  1989).  Area  averages  satisfy 

(_  )  =  (  )„<*  +  (  Ul-cr). 

il) 

where  an  overbar  denotes  an  area  average.  Subscripts  u  and  d 
denote  upward  and  downward  moving  parcels,  respectively. 

Consider  an  arbitrary  scalar  i //.  The  turbulent  flux  of 
i ff  associated  with  the  convective  circulauons  is  given  by 

f„  =  Pvv>'  =  p[(w„-  w)(  +(w3  -  iv)(  t/^  -  H/ )( 1  -  CT  l 


where 

Mtsp<r(l-(j  )(*.-*-,) 

1 3 1 

is  the  convective  mass  flux. 

We  enforce  consistency  between  the  mass  fiux' 
representation  of  the  fluxes  and  the  boundary  conditions  on  the 
fiuxes  at  the  surface  and  the  PBL  top.  For  brevity,  we 
consider  here  only  the  PBL  top.  The  assumption  that  the 
entrainment  layer  is  thin  yields  the  familiar  "jump  relation 
between  <Fv)b  and  the  entrainment  rate: 

=  ¥,)-  js.d:. 

•• 

i  4  i 

Here  we  follow  Lilly  (1968)  by  keeping  the  S,  term,  which 
represents  a  possible  concentrated  entrainment- layer  source 
of  yz,  (e.g.,  due  to  radiation).  Consistency  implies  that 

>> 

My  Vd~V>,=  £<*,<¥,.  -  ¥,) +  cr,  J  S.d: 

(5) 


At  this  point,  we  introduce  a  mixing  parameter  £f.  that 
relates  the  propemes  or'  the  descending  air  at  level  B  to  those 
of  the  free  atmosphere  and  the  mean  state  at  level  B.  We 
include  a  term  that  represents  the  effects  of  the  concentrated 
source: 


~  ¥  >a  =  X e<  Va.  '  ¥3/  +  A  J  S,  dz. 

•« 

(6) 


In  such  dry  adiabatic  cases,  the  entramment-layer  parcels  w  ith 
the  lowest  virtual  dry  static  energies  are  those  for  which  x  is 
zero.  [Cloudy  cases  with  phase  changes  and  radiative  cooling 
are  discussed  later.) 

Let  the  pdf  for  x  in  the  entrainment  layer  be  denoted 

by  nslx).  When  (Hi  is  satisfied,  the  average  virtual  dry 
static  energy  of  the  sinking  air  emerging  from  the  entrainment 
layer  is 


Here  A  is  a  coefficient  that  will  be  determined  later. 


The  mixing  parameter  Xe  1S  closely  related  to  the 
parameter  %  discussed  by  Stems  et  al.  [  1989:  see  also  Albrecht 
et  al.  (1985.1  and  Nicholls  and  Turton  (1986)).  We  can 
interpret  Xe  ns  the  value  of  x  associated  with  the  downdraft  air 
at  level  B.  Since  there  is  a  sharp  gradient  of  y  across  the 
entrainment  layer,  we  expect  0  <  &<<  /. 

Comparing  (3.6)  with  (3.7),  we  find  that 


•  «• 

I-W :,XE  +  VV  +  l-V  SA  *cs)js^d:  =  0. 


(7) 


In  case  the  source  term  of  (3.8)  vanishes,  we  obtain 


fs.(x)nE(x)<*x 

co.dcr,  p«rcen 

J  \x)dx 

„  olden  p*rce_» 


We  assume  that 


l-oa=  \nc(x)dx 

cotden  c**-"ce_j 


L'se  of  (13)  m  (12)  gives 

(04(/_<t*)=  js.(x)nE(x)dx 

coldest  parcel 


(12i 


1 1- 


XeK,  = 


(8) 


This  relationship  does  not  involve  1 \p;  it  must,  therefore,  apply 
for  all  i)/  To  ensure  that  (3  9)  will  be  satisfied  even  when  the 
source  term  of  ( 3.8)  is  not  zero,  we  must  choose 


A  =  a,/  M:l. 


(9) 


We  can  interpret  (8)  as  a  "continuity  equation'"  for  the 
eddies.  We  can  use  (9)  to  eliminate  A  in  (7);  the  result  is: 


¥*  =  Xe  1  -  Ze> Vs  +  j 5, d: . 

(10) 

According  to  ( 10).  the  descending  air  at  level  B  has  the 
propemes  of  the  free  atmosphere,  except  as  modified  by 
mixing  (when  Xe  <  7)  and  by  the  concentrated  source. 
Caughey  et  al.  (1982)  and  Nicholls  (1989)  have  reported 
observations  of  cool  downdrafts  in  the  upper  portions  of 
stratocumulus  cloud  sheets.  They  concluded  that  the  sinking 
air  had  been  radiatively  cooled  near  the  cloud  top.  Such 
effects  are  represented  by  the  5 v  term  of  ( 10). 

The  next  step  is  to  introduce  a  further  simple  relation 
between  X  and  cr.  We  assume  that  the  descending  parcels  in 
the  convective  circulations  at  level  B  are  formed  from  the 
densest  parcels  available  in  the  entrainment  layer.  The  various 
parcels  are  produced  by  mixing,  in  various  proportions,  air 
from  the  intenor  of  the  PBL  with  air  that  has  recently  been 
entrained.  For  the  entrainment  layer,  define  X  without  a 
subscnpt  as  follows:  for  each  entrainment-layer  parcel,  x  IS 
the  mixing  fraction  of  the  air  with  free-atmosphenc  properties. 
Under  dry  adiabatic  mixing,  we  have 


ulA)  =  (/-A)*,,- 


(ID 


Simple  algebra  leads  to 


Ea,(l-at)= 

c*>L0**i  pmej 

*(15) 

Parcels  with  propemes  very  close  to  those  of  the  free 
atmosphere  itself  make  up  only  a  any  fraction  of  all  the  parcels 
available  in  the  entrainment  layer  at  any  given  time.  >0  that 

nE{x)  is  very  small  for  values  of  x  near  one.  and  increases 
significantly  only  as  x  decreases  to  values  near  zero 

Suppose  that  flE{x)  is  zero  for  Xb  -  X  -  1  ■  and  assumes  a 

constant  value  flEiotO<X-X>-  Clearly,  we  must  require 
that 

x,nE  =  i 

(  lb) 

The  PBL  is  normally  capped  by  an  inversion,  so  that 
the  densest  parcels  available  are  those  for  which  x  1S  near 
zero.  A  straightforward  analysis  leads  to 


/ 


(  17, 

If  the  mean  virtual  dry  static  energy  should  happen  to  decrease 
upward  across  the  entrainment  layer,  the  densest  parcels 

available  will  have  x  =  Xt  For  such  a  statically  unstable 
entrainment  layer,  we  find  that 


/ 


A  similar  analysis  can  be  given  for  the  ventilation  laser, 
allowing  for  both  stable  and  unstable  cases.  1  Randall  and 
Shao.  1990). 


Suppose  that  V/.  and  cr  are  independent  of  height.  By 
combining  our  results  for  the  ventilation  and  entrainment 
layers,  we  can  derive  expressions  for  both  Mc  and  <T  .  The 

forms  of  the  resulting  expressions  guarantee  that  \f  >  0  and 

0<a<l  Plots  are  given  in  Fig.  3.  In  all  cases,  a 
decreases  as  r]  increases.  This  means  that  strong  entrainment 
is  associated  with  small  cr  When  the  ventilation  layer  is 
unstable  and  the  entrainment  layer  is  stable,  the  normalized 
mass  flux  decreases  as  17  increases.  The  reverse  is  true  when 
the  ventilation  layer  is  stable  and  the  entrainment  laver  is 
unstable.  When  both  layers  are  unstable,  the  normalized  mass 
flux  is  maximized  when  r\  =  /  When  both  layers  are  stable, 
the  normalized  mass  flux  is  independent  of  r\. 
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Figure  3:  Plots  of  cr  and  the  normalized  mass  flux. 
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UVSE"  denotes  unstable  ventilation  layer  with 
stable  entrainment  layer;  "UVUE"  denotes 
unstable  ventilation  layer  with  unstable 
entrainment  layer;  "SVSE"  denotes  unstable 
ventilation  layer  with  unstable  entrainment  layer; 
and  finally  "SVUE"  denotes  unstable  ventilation 
layer  with  unstable  entrainment  layer. 


where  h  is  the  moist  static  energy,  and  Sh,( / )  represents  the 
effects  of  radiation  on  the  parcel's  moist  stanc  energy.  We  can 

show  that,  if  Oh„(x )  is  independent  of  x  for  the  range  of  x  in 

question,  the  following  simple  relanon  applies  between  oh^X  I 

and  -1 R.  the  area-averaged  radiative  flux  "jump'  across  the 
entrainment  layer: 


=  -  V  gSh' . 

iZ'h 

For  a  given  distribution  of  6 hJx)-  we  can  work  out  syx  I. 
We  then  identify  the  densest  available  parcels,  in  the  usual 
way. 

3.  Comparison  with  large-eddv  model  results 

To  validate  various  aspects  of  the  parametenzar.cn.  we 
have  used  large-eddy  simulations  produced  by  Moeng  :  luxi, 
1986).  We  now  present  some  preliminary  results  of  this 
exercise. 

Since  drizzle  is  not  included  in  the  large-eddy  model, 
total  moisture  is  a  conservative  vanable.  Passive  dye" 
variables  were  earned  as  additional  conservative  tracers.  We 
have  tested  our  assumptions  of  conservative  mixing  in  the 
entrainment  layer  by  using  moisture  and  dye.  First  we  chose 
model  levels  to  correspond  to  "B"  and  B+'\  For  the  cases 
shown  here,  the  levels  used  were  40  and  44,  respectively. 
The  vertical  resolution  used  was  12.5  m.  so  these  levels 
corresponded  to  heights  of  500  m  and  550  m.  By  averaging 
over  the  domain,  we  computed  mean  values  for  the  dye  and 
moisture  variables  at  B  and  B+.  Then,  for  each  individual  gnd 
volume  found  at  levels  41  -  43,  we  computed  the  two  values 
of  X  corresponding  to  the  local  mixing  ratios  of  dve  and 
moisture.  According  to  our  assumptions,  these  two  values  of 
X  should  agree,  and  of  course  all  values  of  x  should  be 
between  zero  and  one.  Fig.  4  shows  that  the  results  are  quite 
satisfactory. 
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Phase  change  processes  in  clouds  can  be  included  in 
the  above  analysis  without  introducing  any  substantive 
changes  other  than  minor  additional  complexity.  When 
radiative  cooling  occurs  in  the  entrainment  layer,  (1 1)  must  be 
modified  to  take  into  account  the  radiatively  induced  moist 
static  energy  decrease  that  each  parcel  experiences,  as  a 
function  of  x  We  van  write 

Hx)  =  x*Z+<i-x r>^ -*-**(*). 
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Figured:  Scatter  diagram  showing  the  values  ot  x 

obtained  from  total  moisture  and  dye  . 
respectively,  based  on  the  large-eddy  simulations 
of  Moeng  ( 1986). 
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Next,  we  used  the  value  of  x  computed  from  the 

moisture  field  to  infer  Sh„\x)  from  (19).  Fig.  5  shows  the 
results,  in  units  of  degrees  twe  have  divided  by  the  specific 

heat  at  constant  pressure).  For  most  values  of  x  ■  ^iJx )  is 
negative,  as  we  would  expect.  The  values  are  near  -0.1  K. 
For  smail  values  or  x  .  however,  we  obtained  positive  values 
of  Oh s{ x )  Further  inspection  of  the  model  results  shows  that, 
although  cooling  dominates,  radiative  heating  does  in  fact 
occur  locally  at  a  significant  number  of  gnd  points  near  cloud 
top.  The  reasons  for  this  are  currently  being  investigated. 
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Figured:  Values  of  d/iff(^)  plotted  as  a  function  of  x  • 

based  on  the  large-eddy  simulations  of  Moeng 
(1986). 


4.  Conclusions 

We  have  presented  what  amounts  to  very  simple 
second-order-closure  model  that  makes  use  of  an  explicit  (but 
highly  idealized)  description  of  a  turbulent  eddy:  the 
convective  circulation. "  In  adopting  this  approach,  we  have 
necessarily  abandoned  any  pretense  that  the  second-order 
closure  is  so  general  that  it  can  apply  to  any  type  of  turbulent 
boundary  layer.  On  the  other  hand,  we  have  benefitted  by 
being  able  to  use  mechanistic  ideas  about  the  dynamics  of 
individual  turbulent  elements,  which  would  be  difficult  or 
impossible  to  express  in  the  purely  statistical  framework  of  a 
conventional  second-order-closure  model. 

At  present,  we  are  analyzing  tethered  balloon  data 
collected  at  San  Nicolas  Island  during  FIRE,  to  determine  the 
applicability  of  the  model  to  the  observations.  We  are  also 
performing  a  further  analysis  of  Moeng  s  large-eddy 
simulations.  Further  results  will  be  presented  at  the 
Conference. 
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